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I describe the current status of the solar neutrino problem, summarizing the ar- 
guments that its resolution will require new particle physics. The phenomenon of 
matter-enhanced neutrino oscillations is reviewed. I consider the implications of 
current experiments - including the SuperKamiokande atmospheric neutrino and 
LSND measurements - and the need for additional constraints from SNO and other 
new detectors. 



1 Introduction 

Part of the interest in neutrino astrophysics has to do with the fascinating 
interplay between nuclear and particle physics issues — e.g., whether neutri- 
nos are massive and undergo flavor oscillations, whether they have detectable 
electromagnetic moments, etc. — and astrophysical phenomena, such as the 
clustering of matter on large scales, the processes responsible for the synthesis 
of nuclei, the mechanism for core-collapse supernovae, and the evolution of our 
sun. This summary addresses one of the oldest problems in neutrino astro- 
physics, the 30-year puzzle of the missing solar neutrinos. This puzzle grew 
out of attempts to test the standard theory of main sequence stellar evolution, 
but has now led to speculations about physics beyond the standard model 
of electrowcak interactions. I will describe the work that defined the solar 
neutrino problem, the likelihood that its resolution is connected with massive 
neutrinos, and the hopes we have for future experiments. 



2 Open Questions in Neutrino Physics 







The existence of the neutrino was first suggested by Wolfgang Pauli in a private 
letter dated December, 1930. The motivation was to solve an apparent problem 
with energy conservation in nuclear [3 decay: the observable particles in the 
final state (the daughter nucleus and emitted electron) carried less energy than 
that released in the nuclear decay. Pauli suggested that an unobserved particle, 
the neutrino, accompanied the decay and accounted for the missing energy. 

Ajiumber of important developments followed Pauli's suggestion. In 1934 
FermiEl suggested a theory of /3 decay that was modeled after electromagnetism, 
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except that there was no analog of the electromagnetic field: the interaction 
occurred at a point. (Apart from the missing aspect of parity violation, this 
was the correct reduction of today's standard model to an effective theory.) 
In 1934 Bethe and Critchfield described the role of (3 decay in thermonuclear 
reaction chains powering the stars 

{A, Z) ^{A,Z-l) + e+ + j/e 

thus predicting that our sun produces an enormous neutrino flux. In 1956 
Cowan and Reines a succeeded in measuring neutrinos emitted by a reactor 
through the reaction 

De + P ^ n + , 

exploiting the positron and neutron coincidence. (The neutron was detected 
by (n, 7) on a Cd neutron poison.) In 1957 the weak force mediating neu- 
trino interactions was found to violate parity maximally. Later experiments 
found that the j/g was replicated twice more in nature - the i'^ and Vr - each 
accompanying a distinct charged lepton, 

Ve <-> V^i,^ A*~ Vt ^ t~ . 

Finally all of this physics was embodied in the standard electroweak model, out 
of which came the prediction of a new neutral interaction mediating neutrino 
scattering. 

Despite all of this progress, a remarkable number of questions remain. We 
now believe neutrinos are massive, but still have no measurement of an abso- 
lute neutrino mass (only mass differences). Many models attribute the puzzle 
of neutrino mass — why these neutrinos are so much lighter than other stan- 
dard model particles — to scales well beyond the standard model, but we lack 
independent experimental tools for probing these scales. We do not now the 
particle-antiparticle conjugation properties of neutrinos: because they carry 
no standard model charges, both the Dirac (distinct antiparticle) and Majo- 
rana (no distinction between particle and antiparticle) possibilities are open. 
An associated question is the existence of nonzero electromagnetic moments: 
magnetic, charge radius, anapole, and electric dipole. No nonzero moment has 
been measured. 

Finally, there are many questions about the role of neutrinos in astro- 
physics and cosmology. We suspect cosmic background neutrinos contribute 
to dark matter and may influence large-scale structure formation. However 
direct experimental attempts to measure background neutrinos have failed by 
many orders of magnitude to reach the expected density. Type II supernovae 
convert approximately 99% of the energy released in the infall into neutrinos 
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of all flavors. Yet only DeS were detected from SN1987A. Supernova model- 
ers predict that neutrinos play an essential role in the explosion mechanism 
and in the associated nucleosynthesis, yet there is disagreement about the suc- 
cess of neutrino-driven explosions. Finally, there is great interest in mounting 
searches for very high energy astrophysical neutrinos that might be associated 
with active galactic nuclei, gamma ray bursts, etc. 

Given all of these open questions, 70 years after Pauli's original suggestion, 
it would be nice to have a few more answers. There is every indication that 
some answers will come with the resolution of the solar neutrino puzzle. 

3 The Standard Solar Modelli 

Solar models trace the evolution of the sun over the past 4.7 billion years of 
main sequence burning, thereby predicting the present-day temperature and 
composition profiles of the solar core that govern neutrino production. Stan- 
dard solar models (SSMs) share four basic assumptions: 

* The sun evolves in hydrostatic equilibrium, maintaining a local balance 
between the gravitational force and the pressure gradient. To describe this 
condition in detail, one must specify the equation of state as a function of 
temperature, density, and composition. 

* Energy is transported by radiation and convection. While the solar envelope 
is convective, radiative transport dominates in the core region where ther- 
monuclear reactions take place. The opacity depends sensitively on the solar 
composition, particularly the abundances of heavier elements. 

* Thermonuclear reaction chains generate solar energy. The standard model 
predicts that over 98% of this energy is produced from the pp chain conversion 
of four protons into ''He (see Fig. 1) 

ip He + 2e+ + (1) 

with proton burning through the CNO cycle contributing the remaining 2%. 
The sun is a large but slow reactor: the core temperature, ^ 1.5 • 10^ K, re- 
sults in typical center-of-mass energies for reacting particles of 10 keV, much 
less than the Coulomb barriers inhibiting charged particle nuclear reactions. 
Thus reaction cross sections are small: in most cases laboratory measurements 
are only possible at higher energies, so that cross section data must be extrap- 
olated to the solar energies of interest. 

* The model is constrained to produce today's solar radius, mass, and lumi- 
nosity. An important assumption of the standard model is that the sun was 
highly convective, and therefore uniform in composition, when it first entered 
the main sequence. It is furthermore assumed that the surface abundances of 
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metals (nuclei with A > 5) were undisturbed by the subsequent evolution, and 

thus provide a record of the initial solar metallicity. The remaining parameter 
is the initial ^He/H ratio, which is adjusted until the model reproduces the 
present solar luminosity in today's sun. The resulting ^He/H mass fraction 
ratio is typically 0.27 ± 0.01, which can be compared to the big-bang value of 
0.23 ± 0.01. Note that the sun was formed from previously processed material. 

p + p — + e^ + z/ p + p + e" + V 
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^ + p ^^He + 7 
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Figure 1; The solar pp chain. 

The model that emerges is an evolving sun. As the core's chemical com- 
position changes, the opacity and core temperature rise, producing a 44% 
luminosity increase since the onset of the main sequence. The temperature 
rise governs the competition between the three cycles of the pp chain: the 
ppl cycle dominates below about 1.6 -10^ K; the ppll cycle between (1.7-2.3) 
TO^K; and the ppIII above 2.4 •lO'^K. The central core temperature of today's 
SSM is about 1.55 -lO^K. 

The competition between the cycles determines the pattern of neutrino 
fluxes. Thus one consequence of the thermal evolution of our sun is that the 

neutrino flux, the most temperature-dependent component, proves to be 
of relatively recent origin: the predicted flux increases exponentially with a 
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doubling period of about 0.9 billion years. 

A final aspect of SSM evolution is the formation of composition gradients 
on nuclear burning timescales. Clearly there is a gradual enrichment of the 
solar core in "^He, the ashes of the pp chain. Another element, '^He, can be 
considered a catalyst for the pp chain, being produced and then consumed, 
and thus eventually reaching some equilibrium abundance. The timescale for 
equilibrium to be established as well as the final equilibrium abundance are 
both sharply decreasing functions of temperature, and therefore increasing 
functions of the distance from the center of the core. Thus a steep ^He density 
gradient is established over time. „ 

The SSM has had some notable successes. From helioseismology la the 
sound speed profile c(r) has been very accurately determined for the outer 
90% of the sun, and is in excellent agreement with the SSM. Such studies 
verify important predictions of the SSM, such as the depth of the convective 
zone. However the SSM is not a complete model in that it does not explain all 
features of solar structure, such as the depletion of surface Li by two orders of 
magnitude. This is usually attributed to convective processes that operated at 
some epoch in our sun's history, dredging Li to a depth where burning takes 
place. 

The principal neutrino-producing reactions of the pp chain and CNO cycle 
are summarized in Table 1. The first six reactions produce f3 decay neutrino 
spectra having allowed shapes with endpoints given by E™^'^. Deviations from 
an allowed spectrum occur for neutrinos because the ^Be final state is 
a broad resonance. The last two reactions produce line sources of electron 
capture neutrinos, with widths ^ 2 keV characteristic of the temperature of 
the solar core. Measurements of the pp, "^Be, and ^B neutrino fiuxes will 
determine the relative contributions of the ppl, ppH, and ppIH cycles to solar 
energy generation. As discussed above, and as later illustrations will show 
more clearly, this competition is governed in large classes of solar models by 
a single parameter, the central temperature T^. Theilux predictions of the 
1998 calculations of Bahcall, Basu, and Pinsonneault □ (BP98) and of Brun, 
Turck-Chieze and Morel B are included in Table 1. 

4 Solar Neutrino Experiments and their Implications 

The first solar neutrino results were announced by Ray Davis Jr. and his 
Brookhaven collaborators in 1968, more than 30 years agoQ. Located deep 
within the Homestake Gold Mine in Lead, South Dakota, the detector consists 
of a 100,000 gallon tank of C2CI4. Solar neutrinos are captured by 

37Cl(z/,e-)3Ur. 
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Table 1: Solar neutrino sources and the flux predictions of the BP98 and Brun/Turck- 
Chieze/Morel SSMs in cm-^g-i. 



Source 


^,nax (MeV) 


BP98 


BTCM98 


P + P 




+ ly 


0.42 


5.94E10 


5.98E10 




.13c + c+ - 


h 


1.20 


6.05E8 


4.66E8 


150 - 


,15-^ + e+ + 1/ 


1.73 


5.32E8 


3.97E8 


17f_ 


+ c+ H 


- u 


1.74 


6.33E6 






^Bc + c+ 4- 


- V 


- 15 


5.15E6 


4.82E6 


^He + p 


^^He + e- 




18.77 


2.10E3 




^Bc + 




- V 


0.86 (90%) 


4.80E9 


4.70E9 








0.38 (10%) 






p + 


+ p ^^H 


+ V 


1.44 


1.39E8 


1.41E8 



As the threshold for this reactfon is 0.814 MeV, the important neutrino sources 
are the ^Be and *B reactions. The ^Be neutrinos excite just the Gamow- Teller 
(GT) transition to the ground state, the strength of which is known from 
the electron capture lifetime of ^'^Ar. The '^B neutrinos can excite all bound 
states in '^'^Ar, including the dominant transition to the isobaric analog state 
residing at an excitation energy of 4.99 McV. The strength of excite-state GT 
transitions can be determined from the (3 decay '^''Ga{(3~^)^''^K, which is the 
isospm mirror reaction to ^'^Cl{v, e~)^'^ Ai. The net result is that, for SSM 
fluxes, 78% of the capture rate should be_due to ^B neutrinos, and 15% to 
"^Be neutrinos. The measured capture rate! 2.56 ±0.16 ± 0.16 SNU (1 SNU = 
10~36 capture/atom/sec) is about one-third the SSM value. 

Similar radiochemical experiments were begun in January, 1990, and May, 
1991, respectively, by the SAGE and GALLEX collaborations using a different 
target, ^^Ga. The special properties of this target include its low threshold 
and an unusually strong transition to the ground state of ^^Ge, leading to 
a large pp neutrino cross section (see Fig. 2). The experimental Jiapture 
rates are 66.6l:?;^+!^:^ (SAGE)i an(L77.5 ± 6.2 ^ SNU (GALLEX) 0. The 
SSM prediction is about 130 SNU EJ. Most important, since the pp flux is 
directly constrained by the solar luminosity in all steady-state models, there 
is a minimum theoretical value for the capture rate of 79 SNU, given standard 
model weak interaction physics. Note there are substantial uncertainties in the 
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^^Ga cross section due to ^Be neutrino capture to two excited states of unknown 
strength. These uncertainties were greatly reduced by direct cahbrations of 
both detectors using ^^Cr neutrino sources. 



500 3/2" 



175 5/2 



1/2" 



71 



Ge 



71 



Ga 



pp Be Cr B 
neutrino sources 



Figure 2: Level scheme for '^Ge showing the excited states that contribute to absorption of 
pp, ''Be, ^^Cr and *B neutrinos. 

Experiments of a different kind, Kamiokande II/III and SuperKamiokande, 
exploit water Cerenkov detectors to view solar neutrinos event-by-event. So- 
lar neutrinos scatter off electrons, with the recoiling electrons producing the 
Cerenkov radiation that is then recorded in surrounding phototubes. Thresh- 
olds are determined by background rates; SuperKamiokande is currently op- 
erating with a trigger at approximately six MeV. The initial experiment, 
Kamiokande II/III, found a flux of neutrinpa of (2.80 ±0.19 ± 0.33) • 
10^/cm^s after about a decade of measurement Ej. Its much larger succes- 
sor SuperKamiokande, with a 22.5 kiloton fiducial volume, yielded thc.result 
(2.42 ±0.04 ±0.06) • lO^cm^s after the first 825 days of measurements B This 
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is about 48% of the SSM flux. This result continues to improve in accuracy. 



These results can be combined to limit the principal solar neutrino fluxes, 
under the assumption that no new particle physics distorts the spectral shape 
of the pp and neutrinos. One flnds 

(/.(^Be) - 

4>{^B) ^ 0A7 4>^^^i^B). (2) 

A reduced ®B neutrino flux can be produced by lowering the central temper- 
ature of the sun somewhat, as ^(^B)'--^ T^^. However, such an adjustment, 
either by varying the parameters of the SSM or by adopting some nonstan- 
dard physics, tends to push the 0(^Be)/0(^B) ratio to higher values rather 
than the low one of eq. (12), 



^B) 



T-'^. (3) 



Thus the observations seem difficult to reconcile with plausible solar model 
variations: one observable {(j){^B)) requires a cooler core while a second, the 
ratio (/)(^Be)/(/)(®B), requires a hotter one. _ 

This physics was nicely illustrated by Castellani et al. lid. These authors 
generated a series of nonstandard models by changing the S-factor for the 
p-|-p reaction, modifying the core metalicity, introducing weakly interacting 
massive particles as a new mechanism for energy transport, etc. The result- 
ing core temperature Tc and neutrino fluxes were then determined, and the 
latter were plotted as a function of the former. The pattern that emerges is 
striking (see Fig. 3): parameter variations producing the same value of Tc pro- 
duce remarkably similar fluxes. Thus Tc provides an excellent one-parameter 
description of standard model perturbations. Figure 3 also illustrates the diffi- 
culty of producing a low ratio of (/)(^Be)/0(^B) when Tc is reduced. This result 
is consistent with our earlier argument and shows that even extreme changes 
in quantities like the metallicity, opacities, or solar age, cannot produce the 
pattern of fluxes deduced from experiment (eq. (2)). 

Is it possible to change the solar model in a way that reduces the ''Be/^B 
neutrino flux ratio? It is appears the answer is no in models where the nuclear 
reactions buriLin equilibrium with plausible cross sections. However Gumming 
and Haxton E£l pointed out that a possible exception was a nonequilibrium 
model in which the solar core is mixed on the timescale of '^He evolution, 
about 10^ years. Thus the pp chain is prevented from reaching equilibrium. 
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Figure 3: The responses of the pp, ^Be, and *B neutrino fluxes to the indicated variations in 
solar model input parameters, displayed as a function of the resulting central temperature 
Tc- From Castellani et al. 



This suggestion has some physical plausibiUty because it aUows the sun to 
burn more efficiently, with a cooler core and enhanced ppl terminations. The 
SSM •^He ppafile is known to be overstable, as was first discussed by Dilke 
and GoughEj. Also the possibility of a persistent convective core powered 
by the '^He gradient has been discussed in the literature. (The SSM core is 
convective for about 10^ years because of out-of-equilibrium burning of the 
CNO Jiiycle.) A strong argument against a mixed core was offered by Bahcall 
et al. EZl, who showed that homogenizing the core of the SSM led to very large 
changes in the helioseismology. While this is a sobering result, this test was 
not done in a self-consistent model. However, there is work in progress to 
test the helioseismology of a more realistic mixed core model - one where the 
^He content of the core, the temperature gradient, the nuclear reaction rates, 
and the luminosity are handled consistently 113. If the helioseismology remains 
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unacceptable, this will rule out the only solar model conjecture for producing a 
reduced ^Be/^B flux ratio, which is necessary to produce fluxes closer to those 
observed. 

However, there is a popular argument showing that no SSM change can 
completely remove the discrepancy with experiment: if one assumes undis- 
torted neutrino spectra, no combination of pp, ''Be, and *B neutrino fluxes fits 
the experimental results wellE£l. In fact, in an unconstrained fit, the required 
^Be flux is unphysical, negative by about 2.5a. This is clearly a strong hint 
that one should look elsewhere for a solution! 

The remaining possibility is new neutrino physics. Suggested particle 
physics solutions of the solar neutrino problem include neutrino oscillations, 
neutrino decay, neutrino magnetic moments, and weakly interacting massive 
particles. Among these, the Mikheyev-Smirnov-Wolfcnstein effect — neutrino 
oscillations enhanced by matter interactions — is widely regarded as perhaps 
the most plausible. 

5 Neutrino Oscillations 

One odd feature of particle physics is that neutrinos, which are not required by 
any symmetry to be massless, nevertheless must be much lighter than any of 
the other known fermions. For instance, the current limit on the Ve mass is < 5 
eV. The standard model requires neutrinos to be massless, but the reasons are 
not fundamental. Dirac mass terms mu, analogous to the mass terms for other 
fermions, cannot be constructed because the model contains no right-handed 
neutrino fields. Neutrinos can also have Majorana mass terms 

vlrriLVL and f^mRi/R, (4) 

where the subscripts L and R denote left- and right-handed projections of 
the neutrino field and the superscript c denotes charge conjugation. The 
first term above is constructed from left-handed fields, but can only arise as 
a nonrenormalizable effective interaction when one is constrained to generate 
ttll with the doublet scalar field of the standard model. The second term is 
absent from the standard model because there are no right-handed neutrino 
fields. 

None of these standard model arguments carries over to the more general, 
unified theories that theorists believe will supplant the standard model. In 
the enlarged multiplets of extended models it is natural to characterize the 
fermions of a single family, e.g., v^, e, u, d, by the same mass scale m^). Small 
neutrino masses are then frequently explained as a result of the Majorana 
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neutrino masses. In the seesaw mechanism, 

M.^fo,"^-y (5) 

Diagonahzation of this matrix produces one Hght neutrino, mjight ^ ''^oi^^), 
and one unobservably heavy, mhcavy ~ i^r- The factor (mjj / mji) is the needed 
smah parameter that accounts for the distinct scale of neutrino masses. The 
masses for the Ve^v^^ and i^r are then related to the squares of the correspond- 
ing quark masses TOu, rric, and mj. Taking mn ~ 10^^ GeV, a typical grand 
unification scale for models built on groups like SO(IO), the seesaw mechanism 
gives the crude relation 

m^^ : m^^ : m^^ ^ 2 • 10"^^ : 2 • 10"'^ : 3 • lO^^eV. (6) 

The fact that solar neutrino experiments can probe small neutrino masses, and 
thus provide insight into possible new mass scales m/j that are far beyond the 
reach of direct accelerator measurements, has been an important theme of the 
field. 

Consider for simplicity just two neutrino flavors. The states of definite 
mass are the states that diagonalize the free Hamiltonian. Similarly the weak 
interaction eigenstates are the states of definite flavor, that is, the accom- 
panies the positron in (3 decay, and the accompanies the muon. There is 
every reason to assume that these two bases are not coincident, but instead 
are related by a nontrivial rotation, 

\i^e) = cos6'„|z^i) -I- sin6'„|i^2) 

= ~sm9^\iyi) +cos6y\iy2) (7) 

where Oy is the (vacuum) mixing angle. |__. 
Consider a I'e produced at time t=0 as a momentum eigenstateE3 

Ht^O)) = \iye) = cos e^liyi) + sin 9^11^2). (8) 
The resulting probability for measuring a i^e downstream then depends on 

= l-sin220,sin2(^^^ ^l-isin2 20, (9) 

where the limit on the right is appropriate for large t. (When one properly de- 
scribes the neutrino state as a wave packet, the large-distance behavior follows 
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from the eventual separation of the mass eigenstates.) If the the oscillation 
length 

ora'^c 

is comparable to or shorter than one astronomical unit, a reduction in the solar 
Ve flux would be expected in terrestrial neutrino oscillations. 

The suggestion that the solar neutrino problem could be explained by 
neutrino oscillations was first made by Pontecorvo in 1958, who pointed out 
the analogy with Kq <-> Kq oscillations. From the point of view of particle 
physics, the sun is a marvelous neutrino source. The neutrinos travel a long 
distance and have low energies (~ 1 MeV), implying a sensitivity to 

W>10"^^ey^ (11) 

In the seesaw mechanism, 5m? ^ rn^, so neutrino masses as low as m2 ~ 10^^ 
eV could be probed. In contrast, terrestrial oscillation experiments with accel- 
erator or reactor neutrinos are typically limited to 5m? > 0.1 eV^. (Planned 
long-baseline experiments, though, will soon push below 0.01 eV^.) 

From the expressions above one expects vacuum oscillations to affect all 
neutrino species equally, if the oscillation length is small compared to an as- 
tronomical unit. This is somewhat in conflict with the data, as we have argued 
that the ''Be neutrino flux is quite suppressed. Furthermore, there is a weak 
theoretical prejudice that 9^ should be small, like the Cabibbo angle. The first 
objection, however, can be circumvented in the case of "just so" oscillations 
where the oscillation length is comparable to one astronomical unit. In this 
case the oscillation probability becomes sharply energy dependent, and one 
can choose 5m? to preferentially suppress one component (e.g., the monochro- 
matic ^Be neutrinos). This scenario has been explored by several groups and 
remains an interesting possibility. However, the requirement of large mixing 
angles remains. 



6 The Mikheyev-Smirnov-Wolfenstein Mechanism 



In order to include matter effects, we first consider vacuum oscillations for the 
more general case 

\u{t = 0)) = ae{t = 0)|i/e) + a^(i = 0)|i^^). (12) 
from which one easily calculates 

f ae \ _ 1 f--5m^cos29y 5m^sm29y\ ( ae 



'd!^ \a^i J ^ 4E \ 5m^ sin29v 5m^ cos 29v )■ ^ -^J 
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We have equated x ^ t, thatis, set c — 1. 

Mikheyev and SmirnovEJ showed in 1985 that the density dependence 
of the neutrino effective mass, a phenomenon first discussed by Wolfenstein 
in 1978, could greatly enhance oscillation probabilities: a I'e is adiabatically 
transformed into a z^^ as it traverses a critical density within the sun. It became 
clear that the sun was not only an excellent neutrino source, but also a natural 
regenerator for cleverly enhancing the effects of flavor mixing. 

While the original work of Mikheyev and Smirnov was numerical, their 
phenomenon was soon understood analytically as a level-crossing problem. 
The vacuum oscillation evolution equation changes in the presence of matter 
to 

.d_/ae \ _ J_ f 2EV2Gfp{x) - Sm"^ cos Sm'^ sin 29v \ f ae 

^di\a^) ^ AE\ Sm^ sin - 2EV2Gfp{x) + 6m'^ cos 2dv ) 

(14) 

where Gf is the weak coupling constant and p{x) the solar electron density. 
The new contribution to the diagonal elements, 2E\/2Gfp{x), represents the 
effective contribution to Af^ that arises from neutrino-electron scattering. The 
indices of refraction of electron and muon neutrinos differ because the former 
scatter by charged and neutral currents, while the latter have only neutral 
current interactions. The difference in the forward scattering amplitudes de- 
termines the density-dependent splitting of the diagonal elements of the new 
matter equation. 

It is helpful to rewrite this equation in a basis consisting of the light and 
heavy local mass eigenstates (i.e., the states that diagonalize the right-hand 
side of the equation) , 

|i/l(x)) = cos9{x)\i^e) - sme{x)\iyf_,) 

Wh{x)) ^ sm6{x)\ve) + cos6{x)\vfj_). (15) 
The local mixing angle is defined by 

sin 29 V 



sm29{x) 



X^{x)+sin^ 29z 



coa29{x) ^ ^=^M= (16) 
X^{x)+sin^ 29v 

where X{x) = 2y/2G p p{x)E / 5m? — 008 26*1;. Thus 9{x) ranges from 9v to 7r/2 
as the density p{x) goes from to oo. 
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If we define 

\iy{x)) = aH{x)\iyHix)) + aL(a;)|z^L(x)), (17) 
the neutrino propagation can be rewritten in terms of the local mass eigenstates 



dx 



an 



A(a;) ia{x) \ f gh 
-ia{x) — A(x) J \ aL 



with the splitting of the local mass eigenstates ck^termined by 

6m^ 



2X{x) 



2E 



X'^{x)+sm^ 2ev 



and with mixing of these eigenstates governed by the density gradient 

E \ V2Gf£p{x) sin 29v 



a{x) 



X2(x) +sin^ 26^ 



(18) 



(19) 



(20) 



The results above are quite interesting: the local mass eigenstates diagonalize 
the matrix if the density is constant. In such a limit, the problem is no more 
complicated than our original vacuum oscillation case, although our mixing 
angle is changed because of the matter effects. But if the density is not con- 
stant, the mass eigenstates in fact evolve as the density changes. This is the 
crux of the MSW effect. Note that the splitting achieves its minimum value, 
sin 2^^, at a critical density pc = p{xc) 



2-^EGfPc = Sm^ cos 20v 



(21) 



that defines the point where the diagonal elements of the original flavor matrix 
cross. 

Our local-mass-eigenstate form of the propagation equation can be trivially 
integrated if the splitting of the diagonal elements is large compared to the off- 
diagonal elements, 



j{x) 



Xix) 



a{x) 



sin^ 2ev Sn? 



1 



{X{xf 



cos26'v 2E I 1 Mx) I 

^ Pc dx I 



sin'^ 2ey 



> 1, (22) 



a condition that becomes particularly stringent near the crossing point, 

sin^ 26^ Sm^ 1 



7c = l{Xc) = 



cos 26*,; 2E 



1 dp(x) I 



» 1. 



(23) 
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The resulting adiabatic electron neutrino survival probability valid when 
7c > 1, is 

Pi^f"""" " ^ + ^ ^^"^^ 

where 6i = 9{xi) is the local mixing angle at the density where the neutrino 
was produced. 

^(x) -7r/2 ^(x) 
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Figure 4: Schematic illustration of the MSW crossing. The dashed lines correspond to the 
electron-electron and muon-muon diagonal elements of the matrix in the flavor basis. 
Their intersection defines the level-crossing density pc. The solid lines are the trajectories 
of the light and heavy local mass eigenstates. If the electron neutrino is produced at high 
density and propagates adiabatically, it will follow the heavy-mass trajectory, emerging from 
the sun as a 

The physical picture behind this derivation is illustrated in Fig. 4. One 
makes the usual assumption that, in vacuum, the is almost identical to the 
light mass eigenstate, i^l(O), i.e., mi < TO2 and cosO^ ^ 1. But as the density 
increases, the matter effects make the Vg heavier than the i/^, with i^e — > vnix) 
as p{x) becomes large. That is, the mixing angle at high density rotates to 
7r/2. The special property of the sun is that it produces v^s at high density 
that then propagate to the vacuum where they are measured. The adiabatic 
approximation tells us that if initially Vf. vh {x) , the neutrino will remain on 
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the heavy mass trajectory provided the density changes slowly. That is, if the 

solar density gradient is sufficiently gentle, the neutrino will emerge from the 
sun as the heavy vacuum eigenstate, ~ u^. This guarantees nearly complete 
conversion of fgS into v^s, producing a flux that cannot be detected by the 
Homestake or SAGE/GALLEX detectors. 

But this does not explain the curious pattern of partial flux suppressions 
coming from the various solar neutrino experiments. The key to this is the 
behavior when 7c < 1- Our expression for 7(.t) shows that the critical region 
for nonadiabatic behavior occurs in a narrow region (for small 9y) surrounding 
the crossing point, and that this behavior is controlled by the derivative of the 
density. This suggests an analytic strategy for handling nonadiabatic crossings: 
one can replace the true solar density by a simpler (integrable!) two-parameter 
form that is constrained to reproduce the true density and its derivative at 
the crossing point Xc- Two convenient choices are the linear {p{x) = a + 
hx) and exponential {p{x) = ae~^'^) proflles. As the density derivative at Xc 
governs the nonadiabatic behavior, this procedure should provide an accurate 
description of the hopping probability between the local mass cigenstates when 
the neutrino traverses the crossing point. The initial and ending points Xi and 
Xf for the artiflcial profile are then chosen so that p{xi) is the density where 
the neutrino was produced in the solar core and p(xf) — (the solar surface), 
as illustrated in in Fig. 5. Since the adiabatic result {P^f^'^^) depends only on 
the local mixing angles at these points, this choice builds in that limit. But 
our original flavor-basis equation can then be integrated exactly for linear and 
exponential profiles, with the results given in terms of parabolic cylinder and 
Whittaker functions, respectively. 

That result can be simplified further by observing that the nonadiabatic 
region is generally confined to a narrow region around Xc, away from the end- 
points Xi and Xf. We can then extend the artificial profile to a; = ±oo, as 
illustrated by the dashed lines in Fig. 5. As the neutrino propagates adiabati- 
cally in the unphysical region x < Xi, the exact soluation in the physical region 
can be recovered by choosing the initial boundary conditions 

ai,(— oo) = — o^(— oo) = cost'^e •'—=^ 

a//(— oo) = ae(— oo) = sinOie' •^-'^^^^^^^ . (25) 

That is, |i^(— oo)) will then adiabatically evolve to \iy{xi)) = |z/e) as x goes from 
— oo to Xi. The unphysical region x > Xf can be handled similarly. 

With some algebra a simple generalization of the adiabatic result emerges 
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Figure 5: The top figure illustrates, for one choice of sin^26 and &m? , that the region 
of nonadiabatic propagation (solid line) is usually confined to a narrow region around the 
crossing point r^- In the lower figure, the solid lines represent the solar density and a 
linear approximation to that density that has the correct initial and final values, as well 
as the correct density and density derivative at Tc- Thus the linear profile is a very good 
approximation to the sun in the vicinity of the crossing point. The MSW equations can 
be solved analytically for this wedge. By extending the wedge to ±00 (dotted lines) and 
assuming adiabatic propagation in these regions of unphysical density, one obtains the simple 
Landau-Zener result discussed in the text. 



that is valid for all Siv? / E and Oy 

Pi,.^\ + \ cos 20^ cos 261,(1 - 2Fhop) (26) 

where Phop is the Landau-Zener probability of hopping from the heavy mass 
trajectory to the light trajectory oil traversing the crossing point. For the 
linear approximation to the dcnsitytJO, 

^"op = e-"^=/'. (27) 

As it must by our construction, P^^ reduces to p^'i'^t) ^^j. ^ 1 When the 
crossing becomes nonadiabatic (e.g., 7c ^ 1 ), the hopping probability goes to 
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1, allowing the neutrino to exit the sun on the light mass trajectory as a Vg, 
i.e., no conversion occurs. 

Thus there are two conditions for strong conversion of solar neutrinos: 
there must be a level crossing (that is, the solar core density must be sufhcient 
to render Vg ~ i^nixi) when it is first produced) and the crossing must be 
adiabatic. The first condition requires that Sni^/E not be too large, and the 
second 7^ > 1. The combination of these two constraints, illustrated in Fig. 

6, defines a triangle of interesting parameters in the ^ sin 20y plane, 

as Alikhcycv and Smirnov first found. A remarkable feature of this triangle 
is that strong conversion can occur for very small mixing angles 

(sin^ 29 10~^), unlike the vacuum case. 




Figure 6: MSW conversion for a neutrino produced at the sun's center. The upper shaded 
region indices thoses 5m? /E where the vacuum mass splitting is too great to be overcome 
by the solar density. Thus no level crossing occurs. The lower shaded region defines the 
region where the level crossing is nonadiabatic (7c less than unity). The unshaded region 
corresponds to adiabatic level crossings where strong z^e — ► will occur. 
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One can envision superimposing on Fig. 6 the spectrum of solar neutri- 
nos, plotted as a function of for some choice of Sm."^. Since Davis sees 
some solar neutrinos, the solutions must correspond to the boundaries of the 
triangle in Fig. 6. The horizontal boundary indicates the maximum for 
which the sun's central density is sufficient to cause a level crossing. If a spec- 
trum properly straddles this boundary, we obtain a result consistent with the 
Homestake experiment in which low energy neutrinos (large 1/E) lie above 
the level-crossing boundary (and thus remain I'e's), but the high-energy neu- 
trinos (small 1/E) fall within the unshaded region where strong conversion 
takes place. Thus such a solution would mimic nonstandard solar models in 
that only the *B neutrino flux would be strongly suppressed. The diagonal 
boundary separates the adiabatic and nonadiabatic regions. If the spectrum 
straddles this boundary, we obtain a second solution in which low energy neu- 
trinos lie within the conversion region, but the high-energy neutrinos (small 
1/E) lie below the conversion region and are characterized by 7 <C 1 at the 
crossing density. (Of course, the boundary is not a sharp one, but is char- 
acterized by the Landau-Zencr exponential). Such a nonadiabatic solution is 
quite distinctive since the flux of pp neutrinos, which is strongly constrained 
in the standard solar model and in any steady-state nonstandard model by 
the solar luminosity, would now be sharply reduced. Finally, one can imag- 
ine "hybrid" solutions where the spectrum straddles both the level-crossing 
(horizontal) boundary and the adiabaticity (diagonal) boundary for small 6, 
thereby reducing the ^Be neutrino flux more than either the pp or fluxes. 

What are the results of a careful search for MSW solutions satisfying 
the Homestake, Kamiokande/SuperKamiokande, and SAGE/GALLEX con- 
straints? One solution, corresponding to a region surrounding Srn^ ~ 6 • 
10~^eV^ and sm^29y ~ 6 • 10~^, is the hybrid case described above. It is 
commonly called the small-angle solution. A second, large-angle solution ex- 
ists, corresponding to Srn^ ~ 10~^eV^ and sin^ 29y ~ 0.6. (Variations on 
these solutions include oscillations to sterile neutrinos, oscillations modified 
by regeneration as neutrinos pass through the earth (day/night effects), and 
of course vacuum oscillations.) These solutions can be distinguished by their 
characteristic distortions of the solar neutrino spectrum. The survival prob- 
abilities P^^^^(E) for the small- and largo-angle parameters given above are 
shown as a function of E in Fig. 7. 

The MSW mechanism provides a natural explanation for the pattern of 
observed solar neutrino fluxes. While it requires profound new physics, both 
massive neutrinos and neutrino mixing are expected in extended models. The 
small-angle solution corresponds to Sm^ ^ 10^^ cV^, and thus is consistent 
with 7712 ~ few -10"^ eV. This is a typical Vt mass in models where niR ~ 
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Figure 7: MSW survival probabilities P(E^) for typical small angle and large angle solutions. 



mcuT- This mass is also reasonably close to atmospheric neutrino values. 
On the other hand, if it is the participating in the oscillation, this gives 
nin ~ 10^^ GeV and predicts a heavy Vr ~ 10 eV. Such a mass is of great 
interest cosmologically as it would have consequences for supernova physics, 
the dark matter problem, and the formation of large-scale structure. 

There are many interesting elaborations of the MSW effect not discussed 
here, but treatedJii many papers: spin-flavor oscillations induced by the solar 
mag netic field BEI (the mass difference between and a sterile lyj^ is com- 
pfipsatfid by the matter effects); oscillations induced by density fluctuations 
EZlCa'E2l; "stochastic depolarization" effects in large random magnetic fields E3; 
etc. There are also interesting effects associated with three neutrinos: if the 



solar neutrino problem is due to MSW Ve 



oscillations, one might expect a 



Vf. — Vt crossing at still higher densities. This has led to many interesting spec- 
ulations about the role of the MSW mechanism in supernova explosions and in 
supernova nucleosynthesis, and about the possibility that v^^ — oscillations 
governed by small mixing angles might be best probed using the supernova 
neutrino flux. 
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7 Other Neutrino Mass Evidence and Implications 

The solar neutrino problem is not the most compelling evidence for neutrino 
mass. SuperKamiokande's analysis of the ratio of muon-like to electron-like 
atmospheric neutrino events confirmed that a dramatic anomaly exists. The 
quantity studied is 

„ _ {N,/N,)data 

[N,/N^)mc ' ^'^^ 

the measured ratio of muon-like to electron-like neutrino events normalized to 
the expected ratio, based on Monte Carlo calculations of the productioii_and 
interaction of atmospheric neutrinos. The SuperKamiokande results areEll 

R ^ 0.63 ± 0.03(stat) ± 0.05(syst) (29) 

for sub-GeV events which were fully contained in the detector and 

R = 0.65 ± 0.05(stat) ± 0.08(syst) (30) 

for fully- and partially-contained multi-GeV events. The results for R are 
consistent among the four largest detectors used for atmospheric neutrinos 
(SuperK, Soudan II, 1MB, Kamiokande). While this suggests neutrino oscilla- 
tions, even stronger evidence for new physics comes from measurements of R as 
a function of the zenith angle, O, between the vertical and neutrino direction. 
A down-going neutrino (O ^ 0°) travels through the atmosphere above the de- 
tector (a distance of about 20 km), whereas an up-going neutrino (0 ~ 180°) 
has traveled through the entire Earth (a distance of about 13000 km). Hence 
a measurement of number of neutrinos as a function of the zenith angle yields 
information about their numbers as a function of the distance traveled. 

The zenith angle dependence of the electron and muon fluxes Lj is shown 
in Fig. 8, plotted as a function of the reconstructed LjE^. The muon neu- 
trino flux drops with increasing distance, while the electron neutrino flux is 
approximately constant. This behavior is consistent with v^j_ <-> i^TpCiscillations. 
Measurements of up-going muons at Kamiokandel23 and MACROc^, produced 
by high energy muon neutrino interactions in the rock below the detectors, 
show a similar deficit. 

It is often pointed out that such results provide convincing evidence be- 
cause the up/down difference in R is essentially self-normalizing: only weak 
geomagnetic effects are expected to break the isotropy of cosmic ray interac- 
tions in the atmosphere. 

Explanations of the SuperK anomaly in terms of ^ Ve oscillations 
conflict with reactor oscillation limits, while the explanation v^^ <-> lygterUe is 
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Figure 8: The ratio of fully contained events measured at SuperKamiokande versus recon- 
structed L/Eu. The dashed lines show the expected shape for j/p <-> v-r oscillations with 
&-m? = 2.2 • 10-3 eV^ and sin^ 29 = 1. 

disfavored in fits. The preferred solution is <-!• v^- characterized by 

5m23 - (4±2) • IQ-^eV'^ 
0.8 < sin^ 26* < 1.0. 

The best fit mixing angle is approximately maximal, an intriguing and surpris- 
ing result. 

One can take the the square root of the atmospheric 5m? to obtain ^ 0.06 
eV, a minimum value for a neutrino mass. This immediately establishes a 
lower bound on the neutrino contribution to dark matter of about 0.3% the 
closure density, a value not too different from the mass evident in visible stars, 
a remarkable result o. 

There is one other indication of neutrino mass, the positive signal for 
— !■ Vf. oscillations seen by the LSND group in a beam-stop experiment 
at Los Alamos Ea. The experiment uses a 52,000 gallon tank of mineral oil 
and liquid scintillator, instrumented with 1220 phototubes. A neutrino event 
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i/e + p — !■ n + e+ is detected by the coincidence of the positron and subsequent 
2.2 McV 7 ray from neutron capture, n + p ^ d + 7. The signal is consistent 
with Df^ — > Pe osciUations in a narrow band that includes the ranges Sm^ 0.2 
- 2.0 eV^ and sin^ 20 ~ 0.03 - 0.003. A similar experiment at the Rutherford 
Laboratory, KARMEN, sees no oscillations, but has lower sensitivity E3. A 
recent combined analysis EZI of the two experiments lowers the confidence level 
of the oscillation claim, but finds a parameter region consistent with both 
experiments. An improved experiment at Fermilab has been approved and 
should yield results in 2002. 

These results have some interesting implications. For example, we dis- 

2 

cussed the quadratic seesaw relation earlier, mjight ~ If we use the atmo- 
spheric Sm^ as a rough guide to the mass (or more correctly the 1/3 mass, 
given the large atmospheric neutrino mixing angle), ~ 0.1 eV, and adopt 
for niD the corresponding third-generation quark mass, mjj ~ rritop ^ 180 
GeV, one obtains mn ~ 0.3 • 10^^ GeV. This is a value reasonable close to 
the supersymmetric grand unified scale of ~ 10^^ GeV, a startling result. It 
has inspired some to hope that current neutrino results are giving us our first 
glimpse of physics at the GUT scale. 

One puzzling aspect of atmospheric, solar, and LSND neutrino results is 
that they require three independent Sm'^s. That is, they do not respect the 
relation 

Sm2i + Snii^ + (5m32 — 0. (31) 

Thus either one of more of the experiments must be attributed to some phe- 
nomenon other than neutrino oscillations, or a fourth neutrino is required. 
That neutrino must be sterile to avoid constraints imposed by the known width 
of the Zq. 

8 Outlook 

The argument that the solar neutrino problem is due to neutrino oscillations 
is, in a sense, circumstantial: this conclusion is derived from combining several 
experiments, no one of which requires new particle physics. There is no direct 
observation of new physics analogous to the zenith angle dependence of the 
SuperKamiokande atmospheric results. For this reason there is great interest 
in a new experiment now taking data in the Creighton nickel mine in Sudbury, 
Ontario, 6800 feet below the surface. The Sudbury Neutrino Observatory 
(SNO) has a central acrylic vessel filled with one kiloton of very pure (99.92%) 
heavy water, surrounded by a shield of 7.5 kilotons of ordinary water. SNO 
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can detect electron neutrinos through the charged current reaction 



(32) 



The Cerenkov Hght from the outgoing electron is then recorded in the array 
of 9800 phototubes that surround SNO's central vessel. The spectrum of pro- 
duced electrons is quite hard, making reconstruction of the energy of the Ug 
easier than in the case of neutrino-electron elastic scattering. Thus the exper- 
imenters may be able to detect distortions of the neutrino spectrum resulting 
from the MSW effect. 

SNO will also study the neutral current reaction 



by detecting the produced neutron either through (n, 7) reactions on salt dis- 
solved in the heavy water or in "^Hc proportional counters. In this way the 
experimenters will obtain an integral measurement of the flux of active neu- 
trinos, independent of flavor. Thus a neutral current signal clearly larger than 
the corresponding signal would show that heavy-flavor neutrinos comprise 
a portion of the solar neutrino flux, providing definite proof of new physics. 

The SNO collaboration is expected to make its first announcement of re- 
sults for the charged current reaction as early as summer, 2000. In future 
years, as the SNO neutral and charged current results become precise and 
as SuperKamiokande continues to amass data, the nature of the solar neu- 
trino problem should become much clearer. The hope is that these results, in 
combination with other solar neutrino results (Borexino, GNO, iodine), with 
new atmospheric and (possibly) supernova neutrino measurements, and with 
precision tests of oscillations at accelerators and reactors, will allow us to com- 
pletely characterize the neutrino mass matrix, providing a window on physics 
well beyond the standard model. 

This work was supported in part by the US Department of Energy. 
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